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ABSTRACT 

We improve earlier Galactic bounds that can be placed on the fraction of dark matter in charged elemental 
particles (CHAMPs). These constraints are of interest for CHAMPs whose mass is too large for them to 
have seen through their electromagnetic interaction with ordinary matter, and whose gyroradius in the galactic 
magnetic field is too small for halo CHAMPs to reach Earth. If unneutralized CHAMPs in that mass range 
are well mixed in the halo, they can at most make up a fraction ^ (3-7) x 10"^ of the mass of the Galactic 
halo. CHAMPs might still be a solution to the cuspy halo problem if they decay to neutral dark matter but a 
fine-tuning is required. We also discuss the case where CHAMPs do not populate a spherical halo. 
Subject headings: dark matter - galaxies: halos - galaxies: kinematics and dynamics - galaxies: magnetic 
fields 



1. INTRODUCTION 

Massive particles with integer electric charge (CHAMPs), 
denoted here by X, were considered as dark matter candidates 
in the late eighties (De Riijula et al. 1990; Dimopoulos et 
al. 1990; Gould et al. 1990; Chivukula et al. 1990). Although 
stable CHAMPs are predicted in some extensions of the stan- 
dard model, astrophysical constraints plus bounds from un- 
derground detectors, from balloon experiments and the lack 
of detection of anomalous hydrogen in the sea water, basi- 
cally rule out CHAMPs as dark matter (see Perl et al. 2001 
and Taoso et al. 2008). In particular, the non-detection of 
heavy water in the sea excludes CHAMPs with masses be- 
tween 10 and 10"^ TeV (Verlerket al. 1992). All the above con- 
straints were derived for the standard flux of particles at Earth 
from the Galactic halo. If magnetic fields prevent the flux of 
CHAMPs to penetrate the Galactic disk, one must reevaluate 
eai-lier bounds (Chuzhoy & Kolb 2008). 

Essentially all X~ should have bound to protons, form- 
ing neutraCHAMPs, which decouple from the photon-baryon 
fluid and drive structure formation prior to recombination. 
NeutraCHAMPs reach Earth unimpeded. Searches for neu- 
traCHAMPs in cosmic rays rule out particles with masses be- 
tween 100 and a few 10^ TeV (Barwick et al. 1990). Nev- 
ertheless, heavy-water searches, cosmic rays searches, and 
constraints from overproduction of ^Li (Berger et al. 2008) 
are only relevant if CHAMPs are singly charged, because for 
other charges, a CHAMP no longer behaves like a proton. 

If a significant fraction of the mass of halos is made up by 
CHAMPs, it may have a strong impact on the observable Uni- 
verse (e.g., Chuzhoy & Kolb 2008). It is therefore important 
to constrain the abundance of CHAMPs in galactic halos. Af- 
ter revising Galactic requirements for CHAMPs to be absent 
in the Galactic disk, we give an upper limit on the abundance 
of CHAMPs in the Galactic halo. 

2. SHIELDING THE DISK WITH MAGNETIC FIELDS 

While neutraCHAMPs have no difficulties to penetrate the 
Galactic disk and the solar wind to reach Earth, the penetra- 
tion of unneutralized CHAMPs may be impeded by the pres- 
ence of magnetic fields. Denote by e the electric charge of 
CHAMPs in units of e, the elementary electron charge. For 



mx > 10 TeV and e < 1, CHAMPs in the Galactic halo be- 
have as a collisionless plasma because the self-collision time 
is > 10^ Gyr and the mean free path is > 10*' kpc. Such a 
plasma consists of charged particles influenced only by grav- 
ity and electromagnetic fields. We will consider first the in- 
teraction of halo CHAMPs with the Galactic magnetic field. 

It is well-known that when charged particles interact with 
a magnetized body, a boundary layer that divides two re- 
gions with different conditions is created (Parks 1991). Thus, 
charged particles will penetrate this boundary by some dis- 
tance before they are turned around by the v'x B force (Figure 
[T]). The boundary layer is formed because of the partial pene- 
tration of the charged particles before they are deflected back. 
The orbits described by negative and positive charged parti- 
cles in the neighborhood of the magnetic boundary are drawn 
in Fig.[T] 

Since the Galactic magnetic field is not perfectly plane- 
parallel and has a nonzero turbulent component, i.e.B = B()+b, 
where Bo is the regular (homogeneous) magnetic field and b 
denotes the turbulent field, the particle motion is determined 
not only by the average magnetic field but also by scattering at 
field fluctuations, a stochastic process which requires the solu- 
tion of transport equations with particle ensembles. Depend- 
ing on the magnitude of these fluctuations, we distinguish be- 
tween weak and strong turbulence which leads to different 
physical phenomena (Kallenrode 1998). Particle propagation 
in turbulent fields can be understood as a diffusive process, 
reason why we consider the spatial diffusion of collisionless 
halo CHAMPs into the galactic disk. The diffusion timescale 
Tdiff accross the galaxy disk thickness H for a halo CHAMP, 
is bracketted in the range: 



(1) 



where D|| and D± are the diffusion coefficients parallel and 
transverse to the mean component of the magnetic field, 
which is observed to be parallel to the disk. Within the 
disk, the magnetic field can be considered static because 
Alfven waves propagate with velocities of the order of the 
Alfven speed ~ 6 km s"', which is smaller than the typ- 
ical velocities of CHAMPs ^ V^vx, where vx — 150 km 
s"' is the one-dimensional velocity dispersion for halo parti- 
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Fig. 1 . — Structure of the magnetic boundary layer formed in the caps of 
galactic disks by the partial penetration of the charged particles before they 
are deflected back. We also illustrate the corresponding orbits of particles in 
the neighborhood of such boundary for both types of charge, that is, negative 
and positive. 

77 = (1 + (Bf) I {b^)T\ and on the rigidity x = 27rrz./Aniax, 
with the Larmor radius defined with respect to the total 
magnetic field and Amax the maximum scale of the turbulence 
- H/2 (Giacalone & Jokipii 1999; Casse et al. 2002). Ob- 
servations of the Galactic polarized synchrotron background 
yield 1 < {b^)/{Bl) < 9 (Fletcher & Shukurov 2001, and 
references therein), implying that 0.5 < < 0.9. Since Tdiff 
scales as the inverse of the diffusion coefficients and those are 
essentially a monotonic function of 77, we use 77 ~ 0.5 in our 
estimate of D± in order to give an upper limit on the diffusion 
timescale. Taking advantage of the numerical result by Casse 
et al. (2002) that D±/{riv) ~ 0.3 for Kolmogorov turbulence 
with 7] = 0.5 and x between 0.05 and 0.4, we find that 
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Therefore, the present configuration and strength of the 
Galactic magnetic field can prevent diffusion of (unacceler- 
ated) CHAMPs across the Galactic disk in the life time of the 
disk for mass particles mx < lO^e TeV. The corresponding 
gyroradius for a mass of 10''e TeV moving at 150 km s"' in a 
field of 5/iG is 0.1 pc. 

It is likely that CHAMPs are accelerated to much higher 
velocities by supernova shocks. The inclusion of this effect 
would give a more constraining upper limit on mx ■ 

3. ENERGY LOSS OF CHAMPS IN THE DISK 

In the preceding section we have seen that halo CHAMPs 
with masses nix < lO^e may have difficulty penetrating the 
magnetized Galactic disk, whereas those inside it stay con- 
fined to the disk. Trapped particles in the disk gain en- 
ergy through electrostatic fields, Fermi acceleration in shock 
waves, and its descendants (e.g., Blandford 1994), and loss 
kinetic energy due to Coulomb scatterings with electrons 
and protons of the diffuse interstellar gas. The dissipation 
timescale is rjis = E/\E\, with E = mv\/2 and 



■In A, 



m^x'x 



(3) 



where is the electron density (« 0.025 cm"^ in the solar 
vicinity) and the Coulomb logarithm has a value of about 
20. CHAMPs may avoid strong cooling if the dissipation 
timescale t^s for CHAMPs trapped on the disk be greater than 
the shock acceleration timescale race- Since race ^ 0.01 Gyr, 
the condition 2tcoo1 > Tacc = 0.01 Gyr implies 



10' 



10* 



Region 1 



Region 2 



Region 3 



0.1000 1.0000 



Fig. 2. — Regions of mass-charge space of astrophysical relevance. Halo 
milli-charged particles in region "1" have no difficulty in penetrating the disk, 
whereas CHAMPs in region "3" would interact strongly with the baryonic 
matter in the disk resulting in strong observable consequences, unless its 
abundance is very tiny. For reference, the excluded region from accelerator 
experiments (AC) after Davidson et al. (2000), is also shown. 

This constraint is valid for any value of e provided that mx 
is larger than the electron mass me ~ 0.5 MeV. We will not 
consider the regime nix < '«« because they are excluded for 
10"'^ < e < 1 (Davidson et al. 2000). In Fig.|2]we plot the re- 
gion of mass-charge space for CHAMPs, combining the con- 
straint derived in ^and Eq. (|4]i. CHAMPs in region "2" can 
be suspended in the halo without penetrating the disk and may 
be very evasive for direct terrestrial detection because the flux 
of CHAMPs reaching Earth may be highly suppressed. If 
CHAMPs are thermal relics, they must also satisfy the uni- 
tary bound mx ^120 TeV (Griest & Kamionkowski 1990). 
Combining this constraint with Eq. (|4]l, we find e < 0.2 for 
thermal relics. 

4. THE GLOBAL MAGNETIC SUPPORT 

For simplicity, let us assume for a moment that the Galactic 
magnetic field is horizontal. As depicted in Fig.[Tl CHAMPs 
with mx < 1 X lO^e TeV execute approximately half a gyro- 
orbit before finding themselves back in the unmagnetized re- 
gion and with velocities directed away from the magnetized 
region. This leads to the plasma being excluded from the mag- 
netized disk. In the boundary layer, a current layer develops 
as a result of a thermal, unmagnetized plasma interacting with 
a magnetized region. It is a classical result that the (kinetic) 
motions of individual particles in collisionless plasmas can 
be reconciled with the role inferred for the pressure in MHD 
(e.g.. King & Newmann 1967; Cravens 1997). For example, 
in a boundary layer like the terrestrial magnetopause, the ther- 
mal (or kinetic) pressure of the solar wind is balanced by the 
pressure of the terrestrial magnetic field and this separates the 
interplanetary magnetic field from the magnetospheric cavity. 

Suppose that CHAMPs cannot penetrate down to z = Z^in 
because of the Galactic magnetic barrier (see Since 
CHAMPs are essentially collisionless, the magnetic field is 
the only agent that may support the weight of the CHAMPs 
in the halo. Integrating the equation of vertical equilibrium 
from z = Zmin to z = oo, assuming zero pressure at z = oo, and 
ignoring the weight of coronal gas, we find 



^'^(Zmin) > j PchKJz, 



(5) 
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mx > 2 X 10 e TeV i50kms"' / where Pb is the magnetic pressure, pchiz) the mass density of 
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CHAMPs and K. the vertical positive gravitational accelera- 
tion. For a spherical dark halo, the contribution of CHAMPs 
to the weight term at the solar vicinity is 



PchK.dz = fpxvl = 1 .7/ X 10"'"dyn cm"^ 



Px 



vx 



^0.01MqPC-V Vl50kms-iy ' 

where / is defined as the mass fraction of dark matter in 
CHAMPs. The weight of the CHAMPs produces a confine- 
ment effect. Interestingly, the observed synchrotron emission 
above the plane in the solar neighbourhood implies that the 
scale height of the magnetic field is greater than what would 
be inferred from the weight distribution of the interstellar 
matter (e.g.. Cox 2005). By requiring that the confinement 
of the magnetic pressure is entirely due to the weight of the 
CHAMPs, an upper value on the abundance of CHAMPs can 
be derived. 

The observed synchrotron emission above the plane in the 
solar neighbourhood indicates that the total magnetic field 
strength is 2-5 /iG at a height of z = 1 kpc (Ferriere 2001 ; Cox 
2005; Gaensler et al. 2008). If we identify Z„i„ as the HWHM 
of the magnetoionic disk ^ 1 kpc (e.g., Kalberla 2003) and 
by equating the magnetic pressure at z = Z^in w 1 kpc with 
the weight term, we obtain the desired constraint on /, once 
adopting the highest magnetic value of 5/iG allowed by ob- 
servations: 

\ -1 

-3 I PX^Q 



f<7x 10" 



Vx 

150 km s- 



(7) 
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This constraint is independent of charge e. This estimate is 
very robust to the precise value adopted for Z^in because the 
magnetic field decays very slowly with z. 

In our derivation, we have assumed that the halo is spher- 
ical. Consider now an oblate isothermal dark halo with axis 
ratio q: 



Px(R,z) = 



AirGaq 



(8) 



where Vc is the asymptotic circular velocity at the equatorial 
plane and a = 7"' arcsin7, with 7 = In this model, 

the velocity dispersion is given, within less than 10%, by vx — 
1.16y^(vc/V2), for flattening 0.05 < ^ < 0.5 (e.g. Gerhard & 
Silk 1996). Even though the velocity dispersion for q<\v& 
smaller than in the spherical case, the weight term changes 
only by ~ 10% as compared to the spherical case, even for 
rather flattened halos {q w 0.5). 

So far, it was assumed that CHAMPs and neutraCHAMPs 
are well mixed. Now, we relax that assumption and search 
for the distributions of neutraCHAMPs and CHAMPs that al- 
low to have the largest mass in CHAMPs. This occurs when 
CHAMPs have an almost zero vertical dispersion and settle 
on to layers of negligible thickness at |z| = Zmin, but neutra- 
CHAMPs populate a spherical halo. The vertical density pro- 
file of CHAMPs is then pch{z) = ^chS(\z\ -Z^^^jl, where E,.;, 
is the surface density of CHAMPs at = Rq. From Eq. (|5]l, we 
obtain that Pg ^ at z = Zmin. Taking Kr~6x 10"'' cm 

s"^ at z = 1 kpc (e.g., Holmberg & Flynn 2004), we find Sc/, < 
6 M© pc"^ or, equivalently, / w T,ch/(2px,QRQ) < 3.5 x 10"^. 
However, it is difficult to justify this additional degree of free- 
dom in the model until a non-gravitational mechanism for so 
efficiently dissipation of CHAMP's energy is satisfactorily es- 
tablished. In addition, it is Ukely that Parker instabilities will 



destroy these cold layers of CHAMPs if they are sustained 
against gravity by magnetic fields. 

Consider now a portion of the disk at larger galactocen- 
tric distances, say R = 2Rq. Following the same procedure 
than in the solar neighbourhood, we need to estimate the to- 
tal magnetic pressure at (2/?Q,Zn,in), which should be respon- 
sible to give support to the halo CHAMPs. The large-scale 
magnetic field may have a scaleheight 5-10 times the scale- 
height of the neutral gas disk, so that we may assume that 
B[)(Zmm) — B[)(z = 0). The random magnetic field is expected 
to be roughly in equipartition with the kinetic energy in the 
turbulence. Therefore, its vertical scaleheight should be sim- 
ilar to that of the gas. If magnetic fields are still a barrier for 
halo CHAMPs, then we may assume that Z^in > H and, con- 
sequently, the magnetic pressure by the random component at 
Zniin is less than 10% the pressure by the random field at z = 0. 
Collecting both contributions, we derive an upper Umit for the 
total magnetic pressure at Zmi„: 

;>.<£^±iif = £l„.o.„,, ,9, 

OTT ovr 

where a = b^/Bfi, with and evaluated at z = 0. The 
ordered magnetic field is difficult to measure in the outer 
Galaxy, but there is evidence that it decays with radius R as 
a power-law between /?"' and R~^, probably as expi-R/Rs) 
with Rb = 8.5 kpc (Heiles 1996; Han et al. 2006). The uni- 
form magnetic field in the solar neighbourhood is 2-ApG, 
depending on the authors (Beck 2002; Han et al. 2006). If 
we generously take a value in the solar circle of 4/iG we 
infer a strength Bq ~ 1.5^G at 2Rq. Assuming a spheri- 
cal dark halo with a mass density at 2Rq of ^ px.0/4, then 
/ pchK,dz = fpx.QVx /4. At 2Rq, our assumption that the halo 
is spherical is a very good approximation (e.g., Belokurov et 
al. 2006; Fellhauer et al. 2006). By imposing pressure balance 
at z = Zjnin (Eq.|5]), the following constraint for / is inferred 

-3/1 ,f,i„^f Px.q/4- A"' / Vx 



/<2x 10"-'(l+0.1a) 



0.0025MqPC- 



150 km s 



(10) 

Other observational estimates assure our generously-taken 
magnetic intense. In fact, data from rotation measurements 
of pulsars suggest uniform magnetic fields of ~ 0.7 pG at 
R = 2Rq (Rand & Lyne 1994), which coincides with the ex- 
trapolation of the fit of radial variation of the regular field by 
Han et al. (2006). 

Beyond 2Rq it is uncertain if supernovae shocks are able to 
clean the disk from CHAMPs. It might be also possible that 
beyond the optical radius, the magnetic field is too weak to 
prevent CHAMPs from crossing the disk, but any more com- 
plicated analysis is useless in the face of such ignorance. We 
conclude that charged particles can be suspended in the halo, 
so that they would be impossible to detect as they never reach 
the Earth. However, the mass fraction of bare CHAMPs in 
the halo must be rather smafl / < (2-7) x 10"l In a X^-X' 
symmetric Universe, neutraCHAMPs may compose a frac- 
tion < (2-7) x 10"^, because the ratio of their relative abun- 
dances may be as much as 10:1 (Dimopoulos et al. 1990). 
Therefore, CHAMP models require three ingredients: neutra- 
CHAMPs, CHAMPS and neutral dark matter. One possibility 
is that / w 0.5 ' at early times and CHAMPs decay to neutral 
dark matter with a lifetime < 2.75 Gyr to reach f ^ 3 x 10"^ 
at present. 

' We adopt / = 0.5 because positive CHAMPs and neutraCHAMPs are 
expected to be approximately in equal numbers. 
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5. DISCUSSION 

5.1. The origin of cores in galaxies 

Our upper limit on / rules out CHAMPs as an explanation 
for the formation of constant density cores in dark matter ha- 
los. The reduction of the central density after they have driven 
the formation of galactic halos would be insignificant. Even 
if all the CHAMPs were depleted from the central parts of the 
galaxies, the rotation velocity in a certain galaxy would suffer 
a negligible change of (0.1 -0.35)% for / ~ (2-7) x 10"^ If 
CHAMPs decay to neutral particles, / was larger in the past. 
In such scenario, the formation of cores may result from the 
evacuation of CHAMPs in the central regions of galaxies by 
supernova shocks, provided that the decay times > 1 .5 Gyr. 
Combined with the constraint that the decay times < 2.75 
Gyr, discussed in our previous section, a fine-tuning of the 
decay time is required. 

5.2. Ram pressure stripping and collisions of galaxy clusters 

Magnetic fields couple CHAMPs with themselves and with 
ordinary matter. This coupling might cause ram pressure 
stripping of both baryonic and dark matter of subhalos and 
satellite systems. Consider, for instance, the collision of two 
galaxy clusters. Estimates for the magnetic field strength 
in clusters range from roughly 1 - 10 /iG at the center and 
0.1-1 /iG at a radius of 1 Mpc, values that correspond to a 
plasma beta for the CHAMPs, /? = STrP,/,/^^ w 2/ x 10^"^, 
a hot plasma. Even in this dynamically weak magnetic field, 
the mean gyroradius for a CHAMP with = lO^e TeV, is 
< 5 pc at the center and < 50 pc at 1 Mpc. The govern- 
ing equations of collisionless hot plasmas were developed 
by Chew et al. (1956), whose theory is known as the Chew- 
Goldberger-Low approximation. This approximation, which 
leads to MHD equations with anistropic pressure, is satisfac- 
tory when the Larmor frequency is large compared to other 
characteristic frequencies of the problem and the mean parti- 
cle gyroradius is short compared to the distance over which all 
the macroscopic quantities change appreciably (e.g., Spitzer 
1962; Schmidt 1966). Therefore, charged massive particles in 
the halo of galaxy clusters can be described in the fluid-like 
anisotropic MHD approximation; in the merger process, they 
would behave as a clump of fluid, experiencing ram pressure 
stripping and drag deceleration similar to the gas component. 
Since CHAMPs should be attached to the gas component, the 



observed offset between the centroid of dark matter and the 
collisional gas of the subcluster in the Bullet Cluster implies 
/ < 1 (e.g., Natai-ajan et al. 2002; Mai'kevitch et al. 2004). 
Although the current lensing data accuracy is not sufficient 
to derive the mass distribution of the subcluster in the Bul- 
let Cluster, the derived mass estimates of the subcluster leave 
little room for dark matter in the gas bullet. 

Galactic halo CHAMPs may also exert ram pressure on the 
gas component of the LMC and its stream due to their contin- 
uous scattering by the intrinsic magnetic field of the LMC and 
the Magellanic stream. For a Milky Way-type halo of ^ 10'^ 
Mq, a fraction / of 3 x 10"^ implies that the mass in CHAMPs 
could be up to 3 X 10^ Mq and the density at 50 kpc of 
0.9 X 10"^ Mq pc"^. Since these values are smaller than those 
required to explain the mass and extension of the Magellanic 
Stream and the size and morphology of the gaseous disk of 
LMC (Mastropietro et al. 2005), we cannot reduce any fur- 
ther our upper limit on / with the current observations of the 
LMC disk and the Magellanic Stream. 

6. CONCLUSIONS 

Whilst the common wisdom holds that dark matter is neu- 
tral and collisionless, it is important to explore the possibility 
of it having nonzero, not necessarily integer, charge. We have 
considered the pressure support of CHAMPs in our Galaxy to 
derive a simple, upper limit on the fraction of CHAMPs and 
milHCHAMPs, / < (2-7) x 10"l If / was roughly constant 
over time, this constraint rules out CHAMPs as the origin of 
the cores in LSB and dwarf galaxies. In the range of astro- 
physical interest, CHAMPs behave like strongly interacting 
(fluid-like) dark matter (SIDM). Thus, they face many of the 
problems attributed to SIDM. As some examples, we have 
discussed the survival of the Magellanic Stream and the mass 
distribution of the Bullet Cluster. Our constraint that the mass 
in CHAMPs in the Galaxy is not larger than the mass of coro- 
nal gas in the halo seems to apply also to galaxy clusters. 
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